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Abstract 

Perturbative NLO and NNLO QCD evolutions of parton distributions are studied, 
in particular in the (very) small-x region, where they are in very good agreement 
with all recent precision measurements of -Ff (-^j Q^)- These predictions turn out to 
be also rather insensitive to the specific choice of the factorization scheme (MS or 
DIS). A characteristic feature of perturbative QCD evolutions is a positive curvature 
of F2 which increases as x decreases. This perturbatively stable prediction provides 
a sensitive test of the range of validity of perturbative QCD. 



1 Introduction 



The curvature of DIS structure functions like F^a;, Q^), i.e., its second derivative with 
respect to the photon's virtuahty at fixed values of x, plays a decisive role in probing 
the range of validity of perturbative QCD evolutions of parton distributions in the small-x 
region. This has been observed recently [1, 2] and it was demonstrated that NLO(MS) 
evolutions imply a positive curvature which increases as x decreases. However, in contrast 
to [1] where this positive curvature was shown to disagree with the data, the conventional 
jull NLO analysis performed in [2] led to the conclusion that no such disagreement prevails. 
It was therefore concluded [2] that the NLO small-x parton evolution equations are not 
challenged by the small-x data on F^. These rather unique predictions provide a check of 
the range of validity of perturbative QCD evolutions. However, the curvature is a rather 
subtle mathematical quantity which a priori may sensitively depend on the theoretical 
(non)perturbative assumptions made for calculating it. The main purpose of the present 
article is to study the dependence and stability of the predicted curvature with respect 
to a different choice of the factorization scheme (DIS versus MS) and to the perturbative 
order of the evolutions by extending the common NLO (2-loop) evolution [2] to the next- 
to- next-to-leading 3- loop order (NNLO) . 

2 Theoretical formalism 

In the common MS factorization scheme the relevant structure function as extracted 
from the DIS ep process can be, up to NNLO, written as [3, 4, 5] 



Fl{x, Q') - F+ (x, Q') + F,,s{x, Q') + F^{x, Q\ ml) 



(1) 



with the non-singlet contribution for three active (light) flavors being given by 
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where a — a{Q'^) = as{Q'^)/4:7r, C2^l{z) — 5(1 — z), C2^^g is the common NLO coefficient 
function (see, for example, [6]) and a convenient expression for the relevant NNLO 2- 
loop Wilson coefficient Cl^^j^ can be found in [3]. The NNLO Q^-evolution of the flavor 
non-singlet combinations = u + u— (d+d) =Uv — dv and = u + u + d + d — 2{s + s) = 
Uv + dy + 4q — 4s, where q = u — d and s = s, is related to the 3- loop splitting function [7] 
P^s'^, besides the usual LO (1-loop) and NLO (2-loop) ones, P^g and P^g'^, respectively 
[3, 8]. Notice that we do not consider sea breaking effects {u ^ d, s ^ s) since the HERA 
data used, and thus our analysis, are not sensitive to such corrections. The flavor singlet 
contribution in (1) reads 



i F2,5(x, Q') = l{ [C^Z + aC^2l + a'Cfl\ ® E + [aC^^j + a^C^^j] ® g] {x, Q') (3) 

with E(x, Q^) = ^q=u,d,s{(l-\-Q) — Uy-\-dv-\-Aq-\-2s, (72^^ = ^J^^)^^ and the additional common 
NLO gluonic coefficient function Cg^j can be again found in [6], for example. Convenient 
expressions for the NNLO cf} and cf^ have been given in [4] and the relevant 3-loop 
splitting functions P^ \ required for the evolution of Q^) and g{x,Q'^), have been 
derived in [9]. We have performed all Q^-evolutions in Mcllin n-moment space and used 
the QCD-PEGASUS program [10] for the NNLO evolutions. In NNLO the strong couphng 
evolves according to da/dhiQ^ = -E|^o/3£a^+2 where Po ^ U - 2//3, Pi = 102 - 38//3 
and P2 = 2857/2 - 5033//18 + 325/754 and the running a{Q^) is appropriately matched 
at Q = mc = 1.4 GeV and Q = mf, = 4.5 GeV. The heavy flavor (charm) contribution 
F2 in (1) is taken as in [2] as given by the fixed-order NLO perturbation theory [11]. 
The small bottom contribution turns out to be negligible for our purposes. Notice that 
a NNLO calculation of heavy quark production is not yet available. For definiteness we 
work in the fixed flavor factorization scheme, given in (l)-(3), rather than in the variable 
(massless quark) scheme since the results for F2 and its curvature remain essentially 
unchanged [2]. 

The choice of a factorization scheme in NLO, other than the MS scheme used thus far, 
might imply similar effects as the additional NNLO contributions in the MS scheme. For 



example, in the deep inelastic scattering (DIS) factorization scheme [12, 5, 6] the Wilson 
coefficients in (2) and (3) are absorbed into the parton distributions, or more precisely into 
their cvohitions, i.e., into the splitting functions. Disregarding for simplicity all NNLO 
contributions, this transformation to the DIS scheme in NLO is achieved via [4, 5] 
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Instead of (2) and (3), the light u, d, s quark contributions to F|' in the NLO (DIS) 
factorization scheme now simply become 



q=u,d,s 



— X 



^ql{x,Q') + lqt{x,Q') 



DIS 



DIS 



+ -xE{x,Q^)j)is + F^ . 



(7) 



The quantitative difference between the NLO (MS) and NLO (DIS) results will turn out 
to be rather small. Therefore we do not consider any further the DIS scheme in NNLO. 

Having obtained the parton distributions q (x, Q^)dis and g{x, Q^)dis from an explicit 
NLO analysis of F2{x,Q'^) in the DIS factorization scheme, one can transform them to 
the MS scheme via (see [13], for example) 



(-) 



1 {x, Q )dis - a 



C^]l^ Vdis + 
.(1) 



Ci\l ® fi-DIS 
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ix,Q^) + Oia' 



where 
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i)+i(9 + 5.) 



{z' + (1 - zf) In ^—^ - 1 + 8z{l - z) 



(8) 
(9) 
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with f — S. This transformation to the MS scheme then allows for a consistent comparison 
of our NLO(DIS) results with the higher-order results obtained in the MS factorization 
scheme. 

3 Quantitative results 

For the present analysis the valence — Uy, dy and sea w — q, g distributions are 
parametrized at an input scale Ql — 1.5 GeV^ as follows: 

xqy{x,Ql) = Ng^'^il-xf-'-il + Cg^y/^+dg^x + egy-^) (12) 
xw{x,Ql) = N^x'^^il-xY^il + Cy.^/x + dy.x) (13) 

and without loss of generality the strange sea is taken to be s = s = 0.5 q. The normal- 
izations Nu^ and N^^ are fixed by Uydx — 2 and dydx — 1, respectively, and Ng is 
fixed via x(E + g)dx — 1. We have somewhat extended the set of DIS data used in [2] 
in order to determine the remaining parameters at larger values of x and of the valence 
distributions. The following data sets have been used: the small-x [14] and large-x [15] 
HI F| data; the fixed target BCDMS data [16] for and using > 20 GeV^ and 
= Q^{^ - 1) + ml > 10 GeV^ cuts, and the proton and deuteron NMC data [17] 
for > 4 GeV^ and > 10 GeV^. This amounts to a total of 740 data points. The 
required overall normalization factor of the data turned out to be 0.98 for BCDMS and 
1.0 for NMC. The resulting parameters of the various fits are summarized in Table 1. 
The relevant small-x predictions are compared with the HI data [14] in Fig. 1, which are 
also consistent with the ZEUS data [18] with partly lower statistics. The present more 
detailed NLO(MS) analysis corresponds to x^/dof = 715.3/720 and the results are com- 
parable to our previous ones [2]. Our new NLO(DIS) and NNLO(3-loop) fits are also very 
similar, corresponding to x^/dof = 714.2/720 and 712.0/720, respectively, although they 
fall slightly below the common NLO(MS) predictions at smaller values of Q'^. It should be 
emphasized that the perturbatively stable QCD predictions are in perfect agreement with 
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all recent high-statistics measurements of the Q^-dependence of F^{x, Q"^) in the (very) 
small-x region. Therefore additional model assumptions concerning further resummations 
of subleading small-a; logarithms (see, for example, [19]) are not required [7, 9]. 

In Figs. 2 and 3 we show our gluon and sea input distributions in (13) and Table 1 
as obtained in our three different fits, as well as their evolved shapes at = 4.5 GeV^ 
in particular in the small-x region. In order to allow for a consistent comparison in the 
MS scheme, our NLO(DIS) results have been transformed to the MS factorization scheme 
using (8) and (9). Note, however, that the gluon distribution in the DIS scheme is very 
similar to the one obtained in NLO(MS) shown in Fig. 2 which holds in particular in the 
small-x region. This agreement becomes even better for increasing values of Q^. This 
agreement is similar for the sea distributions in the small-x region shown in Fig. 3. Only 
for a; ^ 0.1 the NLO(DIS) sea density becomes sizeably smaller than the NLO(MS) one 
shown in Fig. 3. The NLO results are rather similar but distinctively different from the 
NNLO ones in the very small-x region at Q"^ > Q^. In particular the strong increase of 
the gluon distribution xg{x, Q'^) as x ^ at NLO is somewhat tamed by NNLO 3-loop 
effects [9]. 

Turning now to the curvature of we first present in Fig. 4 our results for -Ff (a;, Q^) 
at X = 10"^, together with a global fit MRSTOl NLO result [20], as a function of [1] 



This variable has the advantage that most measurements lie along a straight line [1] as 
indicated by the dotted line in Fig. 4. All our three NLO and NNLO fits give almost 
the same results which are also very similar [2] to the global CTEQ6M NLO fit [21]. In 
contrast to all other fits shown in Fig. 4, only the MRSTOl parametrization results in a 
sizeable curvature for [2] . More exphcitly the curvature can be directly extracted from 




(14) 



F^{x, Q"^) = aQ{x) + ai{x)q + a2{x)q^ . 



(15) 
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The curvature 02 (x) — \dq F^{x, Q"^) is evaluated by fitting this expression to the predic- 
tions for F^x, Q"^) at fixed values of x to a (kinematically) given interval of q. In Fig. 5a 
we present a2{x) which results from experimentally selected g-intervals [1, 2]: 

0.7 < g < 1.4 for 2 x 10"'' <x < 10~^ 

0.7 <q< 1.2 for 5 x 10"^ < a; < 2 x 10"^ . (16) 

It should be noticed that the average value of q decreases with decreasing x due to the 
kinematically more restricted range accessible experimentally. (We deliberately do 
not show the results at the smallest available x — 5 x 10~^ where the g-interval is too 
small, 0.6 < q < 0.8, for fixing a2{x) in (15) uniquely and where moreover present 
measurements are not yet sufficiently accurate [1, 2]). For comparison we also show in 
Fig. 5b the curvature a2{x) for an x-independent fixed g-interval 

0.6 < g < 1.4 (1.5 <Q^ <12 GeV^) . (17) 

Apart from the rather large values of 02 (x) specific [2] for the MRSTOl fit, our NLO and 
NNLO results agree well with the experimental curvatures as calculated and presented in 
[1] using the HI data [14]. Our predictions do not sensitively depend on the factorization 
scheme chosen (MS or DIS) and are, moreover, perturbative stable with the NNLO 3- 
loop results lying typically below the NLO ones, i.e. closer to present data. It should be 
emphasized that the perturbative stable evolutions always result in a positive curvature 
which increases as x decreases. Such unique predictions provide a sensitive test of the 
range of validity of perturbative QCD! This feature is supported by the data shown in 
Fig. 5a. Future analyses of present precision measurements in this very small-x region 
(typically 10~^ ^ x ^ 10~^) should provide additional tests of the theoretical predictions 
concerning the range of validity of perturbative QCD evolutions. 

Finally, the question arises whether the second derivative of F2 with respect to the 
variable g in (15) is indeed dominated by the curvature = ff^F^ / diluQ^Y which is 
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directly related to the evolution equations and to experiment, since d^F^ = d'^F^/dq'^ is 
a linear combination of F^ = dF^/dlnQ'^ = and if = 0{al): 



2 



-K 



+ F^, 



(18) 



with n = 0.5 GeV^/ {Q^ + 0.5 GcV^). In Fig. 6 we show the two contributions in square 
brackets separately taking k — 0.1 which corresponds to choosing = 4.5 GeV^, i.e. 
0' = 1 as an average of our considered fixed ^-interval in (17). The contribution from the 
slope (first derivative) term is indeed strongly suppressed and the curvature is the 
dominant contribution in (18) in the small-x region in NLO as well as in NNLO. Since 
the suppression depends of course on the chosen value for in k wc show in Table 2 
the separate contributions in square brackets in (18) calculated for three typical values of 
g2 in (17) at a fixed value of x = 10"^ in NLO and NNLO. Even at ^ 1.5 QeV^ if 
dominates over reif and therefore (18) represents a rather clean test of the curvature of 
a structure function. 



4 Conclusions 



Perturbative NLO and NNLO QCD evolutions of parton distributions in the (very) small- 
X region are fully compatible with all recent high-statistics measurements of the Q^- 
dependence of if (x, Q^) in that region. The results are perturbatively stable and, fur- 
thermore, are rather insensitive to the factorization scheme chosen (MS or DIS). Therefore 
additional model assumptions concerning further resummations of sublcading small-x log- 
arithms are not required. A characteristic feature of perturbative QCD evolutions is a 
positive curvature a2{x) which increases as x decreases (cf. Fig. 5). This rather unique 
and perturbatively stable prediction plays a decisive role in probing the range of validity 
of perturbative QCD evolutions. Although present data are indicative for such a behav- 
ior, they are statistically insignificant for x < 10~^. Future analyses of present precision 
measurements in the very small-a; region should provide a sensitive test of the range of 
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validity of perturbative QCD and further information concerning the detailed shapes of 
the gluon and sea distributions as well. 
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Table 1: Parameter values of the NLO and NNLO QCD fits with the parameters of the input 
distributions referring to (12) and (13). Here was evaluated by adding in quadrature the 
statistical and systematic errors. 
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Table 2: The separate slope F|' = dF^ / dXnQ"^ and curvature contributions to (18) in 
the MS factorization scheme at x = 10^'^ and for the fixed g-interval in (17) with k = 0.5 
GeV^/(Q^ + 0.5 GeV^). The results are shown for three representative values of Q'^ of this 
interval. Notice that similarly to (15) we used -FKx, Q^) = Ao(x) + Ai(x) In + ^12(3:^) In^ Q^, 
i.e. F^ = Ai + 2A2 In and if = 2A2. 
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10"^ 10"^ 10"^ 10"^ 

X 

Figure 1: Comparison of our various perturbative fits with the HI data [14] at very small- 
X. Our 3-loop NNLO results always refer to the MS factorization scheme. To ease the 
graphical representation, the results and data for the lowest two bins in Q"^ have been 
multiplied by the numbers as indicated. 
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Figure 2: The gluon distributions at the input scale Qq = 1.5 GeV^, corresponding to (13) 
with the parameters given in Table 1, and at = 4.5 GeV^. For a consistent comparison 
with the NNLO and NLO results in the MS factorization scheme, we have transformed 
our NLO-DIS results to the MS scheme using (9) which are denoted by NLO-DIS|^js . 



14 



1.4 



1.2 



1 - 



0.8 - 



0.6 



0.4 - 



0.2 - 







-i — I — I I 1 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I 1 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 



X q 



NNLO 
NLO-MS 
NLO-DIS I MS 



(GeV^) 




_l I I I I I I I I L 



lO"'' 10" 



10"^ 10" 

X 



10"' 10 







Figure 3: As in Fig. 2 but for the sea distribution xq{x, Q^) where q = u — d. The 
NLO-DIS results have been transformed to the MS factorization scheme using (8) which 
are denoted by NLO-DIS|^jg. 
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Figure 4: Predictions for FK^, Q^) at a; = 10"*^ plotted versus q defined in (14). For 
comparison the global fit NLO result of MRSTOl [20] is shown as well. The global 
CTEQ6M NLO fit [21] is very similar to our NLO and NNLO results as can be deduced 
from [2], and the same holds true for the HI fit [15]. Most small-x data lie along the 
straight dotted line [1]. 
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Figure 5: The curvature 02(2;) as defined in (15) for (a) the variable g-intervals in (16) 
and (b) the fixed g-interval in (17). Also shown are the corresponding MRSTOl NLO 
results [20]. The data in (a) are taken from [1]. The NNLO prediction at the lowest 
X- value coincides with the data (full square). 
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X 

Figure 6: The predicted slope F| = ^Ff/^lnQ^ and curvature F2 appearing in (18) 
for the fixed g-interval in (17), with the suppression factor n = 0.1 corresponding to an 
average = 4.5 GeV^ (g = 1). At smallest values of x, the individual upper curves 
always refer to NLO(MS) and the lower ones to NNLO(MS). 
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